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Abstract 

Reduction of the peak heat loads on the plasma facing components is essential for the 
success of the next generation of high fusion power tokamaks such as the International 
Thermonuclear Experimental Reactor (ITER)l. Many present concepts for accomplishing this 
involve the use of atomic processes to transfer the heat from the plasma to the main chamber and 
divertor chamber walls and much of the experimental and theoretical physics research in the fusion 
program is directed toward this issue. The results of these experiments and calculations are the 
result of a complex interplay of many processes. In order to identify the key features of these 
experiments and calculations and the relative role of the primary atomic processes, simple quasi- 
analytic models and the latest atomic physics rate coefficients and cross sections have been used to 
assess the relative roles of central radiation losses through bremsstrahlung, impurity radiation 
losses from the plasma edge, charge exchange and hydrogen radiation losses from the scrape-off 
layer and divertor plasma and impurity radiation losses from the divertor plasma. This anaysis 
indicates that bremsstrahlung from the plasma center and impurity radiation from the plasma edge 
and divertor plasma can each play a significant role in reducing the power to the divertor plates, 
and identifies many of the factors which determine the relative role of each process. For instance, 
for radiation losses in the divertor to be large enough to radiate the power in the divertor for high 
power experiments, a neutral fraction of 10~ 3 to 10~ 2 and an impurity recycling rate of n e T reC ycle of 
~ 10 16 s nr 3 will be required in the divertor. 
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1. Introduction 

If all of the heating power in large fusion experiments such as the International 
Thermonuclear Experimental Reactor (ITER)^ strikes the divertor plates, the peak heat loads on 
plasma facing components will be very high. These loads can be characterised in terms of the heat 

P q 

flux perpendicular to the open flux tubes at the plasma edge, Q, ~ — , v where q™ is 

2A27la 1 + k 2 
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the MHD safety factor, X is the radial decay length of the power at the midplane, a is the plasma 

minor radius, and K is the plasma elongation. For typical ITER parameters of 300 to 600 MW of 

alpha heating, Qy ~ 2000 — 4000 MW/m 2 . The peak heat flux on a divertor plate can be reduced 

B 

by expansion of the flux surfaces in the vicinity of the X point and by the tilt of the field lines — 

Bj 

and inclination of the divertor plate to 20-40 MW/m 2 . Such heat loads are too high to allow the 
development of a sound divertor design because the surface temperatures are too high for the 
thicknesses required for the divertor be able to survive transient events such as plasma disruptions. 
The "wetted area" of the divertor plates, 2 71 R X x 2 (two divertor legs) x 4 (flux surface 
expansion factor) ~ 4 m 2 with a radial decay length, X, for the power of 0.01 cm, is much smaller 
than the total surface area available in the divertor chamber (200 — 400 m 2 ) or in the main plasma 
chamber (1200 m 2 ). The peak heat loads can be reduced to the 0.6 — 4 MW/m 2 range if the power 
can be spread out on the walls of the divertor chamber or radiated from the main plasma. The ITER 
divertor is designed to maximise the role of atomic processes of charge exchange, hydrogen and 

impurity line radiation, ionisation, and elastic collisions between the recycling gas and the plasma 

9 

in the diverted plasma to spread out the heat and momentum^ (Figure 1). 

Conditions where atomic processes have dispersed the heat and momentum have been 
realised on a number of tokamaks, including divertor experiments-^, and limiter experiments ^ , 
but with lower power levels than needed for a next step experiment such as ITER. The challenge is 
to develop a divertor concept where these effects are strong enough to reduce the energy flux on 
the divertor plate by a factor of at least five and preferably to ten or more. There is an active 
program to develop and validate divertor simulations with the data from present experiments and 
use the simulations to analyse and assess divertor concepts for ITER 6 . The modelling results show 
some general trends. In particular, the calculational results indicate that charge exchange losses and 
hydrogen radiation losses are usually relatively small (< 10 — 20%), but that impurity radiation can 
be large if the conditions are optimal. The sophisticated computer models are very complex and 
time-consuming so that it is difficult to obtain a large number of parameter scans when impurity 
radiation is included. Similarly, the results of divertor experiments are complex and the 
experimental run time limited. To identify and evaluate the key issues involved in these 
experiments and calculations, we have assessed the potential role of the major candidate processes 
which might reduce the peak heat loads in the divertor: bremsstrahlung radiation from the central 
plasma, impurity radiation losses from the plasma edge, charge exchange and hydrogen radiation 
losses from the divertor plasma and impurity radiation from the divertor plasma. These processes 
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all rely heavily on atomic collision effects, and we have assembled and assessed the most recent 
atomic data to use in this analysis. A solution of the power exhaust requirement for fusion 
experiments such as ITER will likely involve a mixture of all of these processes, with 
bremsstrahlung from the plasma center and impurity radiation from the plasma edge and the 
divertor playing the major roles. 

The conditions under which impurity radiation in the divertor might be important is 
illustrated in Figure 2. This divertor concept, proposed for ITER, utilises baffles to confine the 
recycling neutral hydrogen and impurity gases in the divertor chamber and facilitate momentum 
exhaust. The louvers are semi-transparent to the neutrals and ensure that the recycling neutral flux 
is relatively uniform on the divertor plasma, and that the fast neutrals formed by charge exchange 
can transfer their momentum to the wall before striking other neutrals. The recycling neutrals are 
ionised in an ionisation front that stretches from the outer baffle to the central dome. The radiation 
zone would be located upstream of the recycling zone and the energy in the divertor is radiated onto 
water cooled louvers parallel to the plasma. The remaining plasma energy falls on dump targets at 
the bottom of the divertor chamber. 

2. Bremsstrahlung Losses 

Bremsstrahlung plays only a small role in present experiments due to high heating power 
densities that present experiments need to balance the transport losses. Experiments on Doublet 
III-D (DIII-D) 7 with 20 MW of heating has a power density of 20 MW / 23 m 3 « 0.85 MW/m^ 
compared to ITER with 300 MW / 2200 m 3 « 0.13 MW/m 3 . The electron temperatures and 
densities are larger and the local transport losses are smaller in ITER so that Bremsstrahlung can 
play a larger role than in present experiments. 

The ratio of Bremsstrahlung losses, P Brem = C B n 2 e Z eff T Y2 , and alpha heating, 



and fee = nB e /n e , and increases with the impurity level and Z. 

For sample operating conditions for ITER with f a ~ 0.05 — 0.15 and fee ~ 0.1 f a , 
0.15 < P Br em / Pa ^ 0.6 with a nominal value of 0.3 for f « 0.1, f Be « 0.01 (Figure 3). 
Bremsstrahlung losses play a minor role in present experiments due to the lower densities and 
temperatures and high heating power densities compared to fusion experiments such as ITER. For 
a 20 MW heating experiment ( ~ 0.8 MW/m 3 ) such as DIII-D with 1% carbon concentration, 
PBrem/Pheat ranges from 1 to 4 % . For ITER with 2% Be (no He) with 50 MW of auxiliary 
heating and no alpha heating, the bremsstrahlung losses would be very comparable to the auxiliary 
heating for temperatures above 7 or 8 keV. 

3. Charge Exchange and Hydrogen Radiation Losses 

Charge exchange is a potential mechanism to transfer the power from the divertor plasma to 




since (<Jv) DT ~ T 2 and decreases 



where f a = n a /n e 
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the side walls. Several factors, however, intervene to limit its effectiveness. Ionisation rates are 
comparable to charge exchange rates except for T e < 3 — 4 eV (Figure 4). For densities of ~ 10 20 
nr 3 (typical of most divertor plasmas), the recycling neutrals are not able to penetrate very far into 
the divertor plasma since, even in the limit where Xi on iz « ^cx> 




1 v v 0.0065 



™>cr<<»U n(10 20 m- 3 ) 



cm 8 for v ~ 10 4 m/s and <ov> -1.5 10" 14 m 3 /s 



(including collisional radiative effects^) is small compared to the dimensions of the plasma. A 
simple model calculation comparing the heat flux lost via charge exchange and the heat flux on the 
plate indicates that charge exchange losses will be relatively small for realistic parameters from the 
divertor plasma between the X point and the divertor plate. The heat flux on the plate Q p i a te 

(MW/m) per unit toroidal length (Figure 5) is approximately Q plate = n d v d (T d y + 20eV)A 

B tor 

where nd, Td and Vd are the plasma density, temperature and velocity at the divertor plate, 20 eV is 
the sum of the hydrogen ionisation potential with some allowance for hydrogen radiation, A is the 
poloidal width of the divertor plasma and B po i and B tor are the poloidal and toroidal magnetic fields 
in the divertor. The charge exchange losses per unit toroidal length can be estimated as 



'eff j 



where n G and v G are the neutral density and 



n v 3 (21 
Q cx = -y-2l-r H r G where G = min 

3 

velocity, ^ T u is an estimate of the upstream average ion kinetic energy (which could be larger than 
3 

2 T u if the upstream flow speed is near sonic or super-sonic), / is poloidal height of the divertor 

plasma and r is the plasma "reflectivity", the fraction of the neutrals incident on the plasma which 
come back due to charge exchange. Penetration effects are included by estimating the fraction of 
the divertor width which the neutrals can penetrate and weighting the losses by the factor G. The 

reflectivity can be estimated as r ~ - — - — ^^ cx . n v and ndVd are linked by particle flux 

{Ov) rY +(<Tv). . ,. 

\ I CX \ I ionization 

ti v B 

balance: -^(1 - r)2l= n d v d —A where (1-r) is the fraction of particles ionized and r is the 
4 B t 

fraction reflected by charge exchange collisions. Using particle flux balance, the fraction of energy 

Q — T Gv 

lost to charge exchange is f ry = — = . 2 " r . Aeff can be evaluated 

Qcx + Q^e \T u Gr + (T d y + 20 eV)(l-r) 



assuming pressure balance along the field lines: n s T s = n u T u => X eff ~ 



T„ v„ v 



o o 



This analysis (Figure 4) indicates that the fraction of energy lost by charge exchange is 
about 15% or less for T u < 180 eV for ITER conditions (n s = 10 20 nr 3 , T s = 260 eV, A = 0.1 m), 
and about 10 % for the most extreme conditions where T u ~ T s , with similar results for DIII-D ( n s 
= 5 x 10 19 nr 3 , T s = 88 eV, A = 0.1 m). The chief factor reducing the charge exchange losses is 
the poor penetration of the neutrals to regions with high ion temperatures. Pressure balance leads to 
an increase in the density as the temperature falls toward the divertor (except near the divertor 
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plate). The resulting increase in density increases the ionisation rate, and decreases the ionisation 
and charge exchange mean free paths. 

Hydrogen radiation is also a potential energy loss from the scrape-off layer between the 
X point and the divertor plate. However, several effects limit the amount of hydrogen radiation. 
Unlike impurities, only recycling neutral hydrogen atoms from the wall will radiate and they must 
penetrate into the SOL plasma. At high densities (> 10 20 nr 3 ), hydrogen ionisation is enhanced 
and hydrogen radiation is suppressed due to multiple collisions with electrons^. The relative 
magnitude of hydrogen radiation losses due to the neutral flux on the sides of the divertor plasma 
(Fig. 2) can be estimated using our charge exchange loss model by replacing 3/2 T u , the maximum 
energy lost per charge exchange event, by the radiated loss per ionisation. 

f„ rnd = ^ = ^ . For the reference ITER and DIII-D cases (n s = 10 20 

J H-rad + E fi + Tj + U .6 eV 

nr 3 , T s = 120 and 260 eV, and A ~ 0.1 m, respectively), the lack of penetration and the lower 
emissivity at higher densities limit the hydrogen radiation to a few percent (Figure 6). The poor 
penetration of the neutral hydrogen and the suppression of the hydrogen radiation each play a role. 
An additional limit on the effectiveness of hydrogen radiation is imposed by the finite opacity of the 
neutral cloud for the dominant lines of hydrogen which reduces the flux of hydrogen radiation as 
well as increases the effective ionisation rate, Absorption lS a/n o (10 20 nr 3 ) where a=0.002 m for 
L a and .0004 m for H a . A simulation of model DIII-D conditions indicates that for densities in the 
10 20 nr 3 range, the flux of hydrogen radiation can be reduced by a factor of 2 or more 

4. Impurity Radiation from the main plasma edge 

Impurity radiation from the main plasma and from the divertor plasma has the potential to 
spread out the heating power over the main chamber and divertor chamber walls and thereby 
reduce the peak heat loads. The temperature and density ranges for these impurities are in the 1 eV 
to 5000 eV and 10 19 to 10 22 nr 3 range. The progress in atomic physics during the last 15 years 
now allows relatively accurate calculation of impurity ionisation, recombination and excitation rate 
coefficients including direct and indirect ionisation and very detailed treatments of dielectronic 
recombination and excitation using collisional-radiative models which can treat meta- stable levels 
and include density effects ^. Detailed calculations of the impurity emission rates for Be, B, C, Ne 
and Ar including collisional radiative effects have recently been carried out . These new rates are 
more accurate for the nearly neutral species that exist at temperatures of several 100 eV than the 
rates that have been previously used before from the ADPAK code^>14 (Figure 7). The ADPAK 
rates are based on energy levels derived from a screened hydrogenic model and employ scaled 
oscillator strengths and recombination rates so the differences are to be expected. The differences 
are even larger at low temperatures for higher Z elements such as Krypton. To assess the potential 
role of Krypton as an impurity feed gas, we have extended our previous calculations of the lower Z 
elements^ to Krypton for T e from 1 to 200 eV. We find that the ADPAK rates for Krypton below 
100 eV are up to 100 times larger than the more accurate calculations, but that the ADPAK rates 
above -100 eV are reasonably accurate because Kr is then sufficiently ionized that the hydrogenic 
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models used in ADPAK are adequate. We have therefore used our detailed calculations for 
Krypton up to 200 eV and the ADPAK rates for T e greater than 200 eV (Figure 8). 

Many present experiments are able to radiate almost all of the heating power from the main 
plasma edge inside the last closed flux surface^. While edge radiation can be important in 
reducing the energy into the divertor, there are some potential drawbacks to exhausting all of the 
power by edge radiation. Large radiation losses from the main plasma will increase the heat flux on 
the first wall which is already near the engineering limits (~ 0.5 MW/m 2 ) for components that can 
only be replaced infrequently. If the edge density is not sufficiently high, the plasma volume 
required to radiate the power will be large, potentially requiring an increase in the minor radius. 
Large amounts of edge radiation may also adversely affect confinement by reducing the power 
across the separatrix below the threshold needed to reach the H-mode. Using the scalings 
developed from the ITER H-mode threshold database, 

P H _ threshold = 0.025n° 15 B T S or 0An e B T R 25 where n e is the line averaged density, Bx is the 
toroidal field, S is the plasma surface area, and R is the major radius in units of MW, 10 20 nr 3 , T, 
m 2 and m^, the power needed across the separatrix is 100 to 400 MW for ITER 
conditions(depending on the edge conditions), which implies that at least 100 MW will need to be 
exhausted by the divertor. Also, a cool edge may lead to low densities in the divertor, thereby 
increasing the pumping requirements for He exhaust. 

The issues can be characterised using a simple model for radial energy transport at the 
plasma edge (r ~ a) 1^-1 8 : 

Q ± = K^J^ = -nn L(T)=>Q^ = -n 2 fL(T)K^ 

or or or ^ (1) 



Ql~2[ c n 2 J z L z (T e )KdT e 



To simplify the analysis and obtain a range for the achievable values of Qj_, we can assume that 
n e (r) is roughly constant near the edge, or at least has an average value n e , and that k ~ constant 

which is roughly consistent with present experiments (e.g.^ ). This scaling is similar to the 
"INTOR" scaling since n%^ ~ k so that % ~ n _1 . This is roughly consistent with the observation 

that % peaks at the plasma edge as n e drops. With H-mode operation in DIII-D in the scrape-off 

layer near the edge, % ~ 0.25 — 0.5 m 2 /s ' , while for the L-mode or ELMy H-mode conditions, 

% can be as large as 2 m 2 /s(c.f. ^0) For edge densities between 5 and 10 x 10 20 nr 3 ,and % 

between 0.25 and 2 , ~ 0.125 — 2 x 10 20 nr 1 s _1 . With these assumptions, the integral 

2 2 C Tc 

becomes: Q ± ~ 2 n e fjc L 7 (T e ) dT e , where the upper limit on the temperature integral is the 

JO 

temperature on the inside boundary of the radiating layer, which for ITER would be in the 2 to 5 
keV range. This scaling has Qj_ ~ n e V(Kf z ), weaker than the volume loss rate with P ~ n e 2 f z . The 
integral has been evaluated for Be, B, C, Ne and Ar using a collisional radiative model, for Fe 
using ADPAK data and for Kr using a mixture of a detailed atomic model and ADPAK data 
(Figure 9) The integral ceases to increase with T e as each impurity becomes fully ionized, so 
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that higher Z impurities such as Fe and Kr radiate more strongly than lower Z impurities. The 
advantage of the higher radiation rates, however, is off-set by the lower allowed fraction of the 
impurity. 

The surface area of the ITER plasma is ~ 1200 m 2 so that for ITER is 0.25 — 0.5 

MW/m 2 for heating powers of 300 to 600 MW reaching the plasma edge. Assuming that the edge 
density is 5 — 10 x 10 19 nr 3 and that the impurity fraction is 1/3 of the "fatal fraction" for which 
the impurity radiation losses equal the alpha heating power , candidate low Z impurities radiate 
only about -0.1 MW/m 2 (Table 1) and medium Z impurities up to ~ 0.3 MW/m 2 for sufficiently 
high edge densities and high edge conductivities. 

Table 1 Radiation from the plasma edges ( Qj_ ~ 0.2 to 0.4 MW/m 2 required for ITER) 





Be 


C 


Ne 


Ar 


Fe 


Kr 


fatal fraction 


0.14 


0.07 


0.025 


0.0054 


.0027 


.0017 


1/3 fatal fraction 


0.05 


0.023 


0.008 


0.002 


.0009 


.0006 


Q ± /(He (K ± f z )0-5)MWm* 


0.1 


0.4 


1.0 


3.0 


7 


10 


Q (MW/m 2 ) 

for 5 x 10 19 < n e < 10 20 
and 0.125 < K<2 


0.004 to 
0.03 


0.01 to 
0.085 


0.016 to 
0.13 


0.02 to 
0.19 


0.04 to 
0.3 


0.05 to 
0.35 



*K in 10 20 m- 1 s- 1 , n e in 10 20 m-3 

On TEXTOR, Q^'s of 0.1 MW/m 2 have been radiated with stable condition and little or no 

confinement degradation^. Our analysis assumes that the impurities are in coronal equilibrium. 
Rapid recycling of impurities from the limiter and wall in TEXTOR plays a large role in enhancing 
the radiation over coronal equilibrium values. The ITER divertor is designed to localise the 
recycling of neutrals and impurities in the divertor chamber, so that one cannot be assured that the 
same level of enhancement will occur in ITER. Edge impurity radiation will thus play an important 
role in transferring some of the heating power to the first wall but, for the reasons outlined above, 
at least 100 MW or more will cross the separatrix to the scrape-off layer and divertor. 

5. Impurity Radiation from the Divertor Plasma 

The ITER divertor concept is based on the use of impurity radiation in the divertor chamber 
to transfer the power from the divertor plasma to the divertor chamber walls (Figure 2). Such 
"detached" operation has been produced on many experiments divertor experiments by a 
combination of gas puffing and injection of gaseous impurities such as Ne or Ar*. Detached 
conditions appear to be brought about by impurity radiation from the edge plasma which lowers the 
temperature and increases the density in a condensation instability to form a "MARFE-like" state, 
followed by momentum loss by charge exchange and elastic collisions which allow neutral atoms 
to transfer the plasma momentum to the walls(c.f. ^). The plasma temperature must be reduced to 
~ 5 eV for charge exchange and elastic collisions to be significant. 
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The conditions needed to radiate the required energy in the divertor can be assessed in a 
similar fashion to edge radiation by changing the form of the thermal conductivity to account for 
parallel heat conduction and using pressure balance along the field lines. The equations become: 

^ = -n e n z L z (T e ) M Q, = -K T e 2 ' 5 ^ M p e = n e T e • 



2 dx 



° K- f 1 T 2 5 ^ T e _ „2,^ ^ r ^0.5 7, 
2 K oJz^Z 1 e 



1 



d&~ P ; s K f z L z T e °- 5 dT e 



dx 
AG, 



r e oJ z Z e 



where K„ ~ 



3.1xl0 9 ' 



Z eJf \nA 



erg 



cm s eV 



3.5 



andF z (f z ) = 



_ / z (%) 



L z (T e )T e 05 dT e 



f z (%) 



^eff 



1 + 0.01/ Z (%)Z(Z-1) 



wzY/z /c„ = K o Z eff , n e {cm 3 ), T e (eV), L z ( ergs cm 3 * : ) and A<2 y 



ergs 



5 cm 



(2) 



AQ ( ) | — if 

Inpracticalunits : hf J" id j= « 2.5 x 10 5 JtH " L z (T e )T e 05 dT e 
n es \10 m J^JF z V Jo 

where the subscript s denotes the value on the separatrix. This integral has been evaluated for six 
impurities, Be, B, C, Ne, Ar and Kr (Figure 10). 

Two candidate machines have been analysed, a 20 MW DIII-D case and a 1.5 GW fusion 
power ITER case (Table 2). T s has been determined from parallel heat conduction (equation 3). 

T div X5 = — rq n dx « 1 lA with L « \Anq w R 



2.5 



dT 



>x 



dT 



3.5 



1 K ° dx 



-MT 



3.5 



T es « 82.2 eV (q n (GW/m 2 ) R(m) q v Z eff f Mwhere K o « 



3.1xl0 y 



(3) 



Z _ In A 

eff 



erg 



cm s eV 



3.5 



If the radiation is to be entirely located in the divertor, the upper temperature to be used in the 
radiation integral (Eq. 2) is the temperature at the X point, for which L x ~ 0.3 L tot . 

Eq. 2 has a different structure than Eq. 1 for edge radiation. AQn is proportional to T s so 
that larger machines with larger connection lengths have greater radiative capability, but the fusion 
power in larger machines is also larger. Equations 1 and 2 describe the limits on the effectiveness 
of exhausting the power by radiation imposed by heat conduction. The volume radiation loss rate 
scales as n e 2 f z and Q P/R n e Vf z , a much weaker dependence. Eq. 1 has Q [f z ] - 5 , whereas 
Eq. 2 has Q [f z /Z e ff] - 5 which further reduces the radiation efficiency. In addition, the integral in 
Eq. 2 is weighted by T - 5 . 
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Table 2 Typical edge temperatures, connection lengths, and parallel heat fluxes 

P a (MW) Qn(GW/m 2 ) A^m 2 ) L(m) R (m) T s (eV) T x ( eV ) 

DIII-D 20 MW 20 0.47 0.043 22 1.67 120 85 

ITER 1.5 GW 240 1.5 0.16 100 8.00 260 185 

T s may be higher if the heat flow along the field lines is "flux-limited" . The temperature at which 
"flux limiting" become important can be determined by estimating the flux limited heat flow as 

fluta-iim* = W e v e T e « 6.7 x lO- 4 ^^)?; 1 5 (eV) for y ~ 0.1 23 . By equating this to the 

conducted flux in Eq. 3, the condition for the heat flow to be flux limited is approximately 

T e > 58 eV\ " e 3 RqZ eff . For the DIII-D and ITER conditions, T flux _ii mi t ~ 185 and 400eV 
VlO m J 

respectively, well above the temperatures calculated in Table 2 using heat conduction. Thus the 
heat flux is not flux limited, especially near the divertor. 

The divertor radiation capability of six impurities Be, B, C, Ne, Ar and Kr (Figure 10) 
indicates that coronal equilibrium rates are inadequate to radiate 20 MW for DIII-D conditions and 
240 MW for ITER conditions from the entire scrape-off layer or in the divertor below the 
X-point(Table 3). An important assumption is that the impurity concentration is uniform 
throughout the plasma. Strong plasma flows in the divertor plasma would tend to retain and 
compress the impurities in the divertor chamber. However, the present picture in which the 
radiating region is upstream of the recycling region would lead to very low plasma flow velocities 
there so that the thermal force would tend to force impurities toward the main plasma. In addition, 
plasma flows due to drifts and other effects, and turbulence due to ELM's would tend to lead to 
mixing of the impurities. 

Table 3 Comparison of Be, C, Ne, Ar and Kr divertor coronal equilibrium radiation efficiencies 
for DIII-D and ITER. 



Element 


Be 


C 


Ne 


Ar 


Kr 


0.33 x fatal f z (%) 


4.7 


2.23 


0.8 


0.18 


0.09 


V(f z (%)/Z eff ) 


1.7 


1.16 


0.68 


0.37 


0.27 


QlDm-D/V(fz(%)/Zeff) at 120 eV 


0.03 


0.09 


0.3 


0.5 


0.8 


QllDIII-D for 120 eV and n e = 10 20 nr 3 


0.052 


0.10 


0.20 


0.19 


0.22 


QllDm-D/V(f z (%)/Zeff)at85eV 


0.018 


0.05 


0.2 


0.25 


0.3 


QllDIII-D for 85 eV and n e ~ 10 20 nr 3 


0.031 


0.058 


0.14 


0.094 


0.081 


QllITER/V(fz(%)/Z e ff) at 260 eV,n e -10 2 m -3 


0.08 


0.2 


0.7 


2 


4 


QllITER (coronal equilibrium) 


0.14 


0.23 


0.48 


0.75 


1.08 


QllITER/V(fz(%)/Z e ff) at 185 eV,n e ~10 2 °m-3 


0.04 


0.12 


0.4 


0.8 


1.4 


QllITER (coronal equilibrium) 


0.07 


0.14 


0.27 


0.3 


0.38 
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The radiation level is proportional to ^ JJ ' L z T 5 dT so that enhancement of the emissivity 

L would increase the radiative losses. L can be enhanced by charge exchange recombination and 

impurity recycling by altering the ionisation balance toward lower charge states which radiate more 

strongly 16,24-26 jfo e charge exchange recombination rate coefficient, <Gv>cxn ni = 

<ov> cx n e ni (n G /n e ), scales like electron-ion recombination and ionisation with an extra factor of 

no/n e , the fraction of neutral atoms to the electron density, so that the charge state distribution, and 

therefore L(T), can be parametrized by no/n e . Similarly, because the ionisation and recombination 

rates have the form <ov>n e n z 1+ , n e can be factored out to produce an equation for the ionisation 

balance of the form 3n z 1+ /3(n e t)=F(<ov>n z 1+ ), so that the effect of rapid recycling can be 

characterised by n e T reC ycle- The magnitude of effect can be estimated by integrating this equation 

with the initial conditions: n z 0+ (n e t=0)=n tot , n z i+ =0 for i=l,Z, from t=0 to t=T reC y C ie> an d 

parameterizing the radiation equilibrium and emissivity as a function of T e for different values of 
97 

n e Trecycle ■ Using the same formalism as with coronal equilibrium, one can evaluate the n G /n e 
and n e T reC ycle required to enhance the radiation rate to radiate the heating power. Approximately 
n G /n e ~ 10" 2 to 10~ 3 and n e T reC ycle ~ 5 x 10 15 to 5 x 10 16 nr 3 s is required for impurities (with n s ~ 
10 20 nr 3 ) to the power in the scrape-off layer for both the DIII-D 20 MW case and the ITER 240 
MW case (Figure 1 1 — 12, Table 4). Neon appears to be the optimum impurity in that it has the 
least demanding requirements, i.e. the lowest no/n e and the largest n e T reC ycle- 

Table 4 Comparison of the requirements for enhancement of Be, C, Ne, and Ar divertor radiation 
efficiencies for DIII-D and ITER. 



Element 


Be C Ne Ar 


no/n e required to radiate 20 MW (DIII-D) 
n e x recvcle ( s nr 3 )required to radiate 20 MW (DIII-D) 


5 x lO" 2 lO" 2 lO" 2 5 xlO" 2 
10 16 10 16 3xl0 16 5xl0 15 


njn & required to radiate 240 MW (ITER) 
n e x recvc i e (s nr 3 ) required to radiate 240 MW(ITER) 


7 x 10" 3 8 x 10" 3 10" 3 4 x 10" 2 
1016 1Q16 4xl0 16 6xl0 15 



6. Summary 



Transferring the energy from the plasma to the plasma facing components by atomic 
processes is a promising approach for reducing the peak loads for the next generation of fusion 
experiments such as ITER. To better understand the results of the complex simulations and 
experiments, we have used simple models to examine the potential role of central radiation losses 
through bremsstrahlung, radiation losses from the plasma edge, charge exchange hydrogen 
radiation and losses from the scrape-off layer and divertor plasma and impurity radiation losses 
from the divertor plasma. The simple analysis gives results which are consistent with the main 
featurea of the complex models. The conclusion of the simple models is that each process can 
contribute to the solution of the problem, with major roles being played by Bremsstrahlung from 
the plasma core and impurity radiation from the plasma edge and divertor plasma with smaller 
contributions due to charge exchange and hydrogen radiation losses. To achieve the required levels 
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of impurity radiation in the divertor will likely require enhancements of the impurity radiation due 
to charge exchange recombination with no/n e ~ 10~ 2 to 10~ 3 and impurity recycling levels of the 
order of n e T recyc i e ~ 5 x 10 15 to 5 x 10 16 nr 3 s. Neon appears to be the optimum impurity in that it 
requires the lowest n G /n e and the largest n e T reC ycle- 
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Figures: 
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Figure 1. Schematic illustration of energy loss mechanisms in tokamaks due to atomic processes. 
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Figure 2. Schmatic Illustration of Dynamic Gas Target Divertor. 
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P„ /P for ITER as a function 
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Figure 3. PBremsstrahlung / Pa for ITER as a function of Temperature for various values of fne 
with fee = 0. 1 fne-The nominal operating range for ITER (8 keV < <T> < 12 keV, and 
0.05 < fn e < 0.15 is shaded. Value of PBremsstrahlung / Pheat are shown for 20 MW of heating in 
DIII-D (f c = 0.01) and 50 MW of heating in ITER (f Be « 0.02). 
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Figure 4. The fraction of charge exchange losses f cx for DIII-D and ITER conditions, <av> for 
ionization and charge-exchange for n e =10 20 nr 3 and <gv>cx/(«3v>cx- +«^v>i on ization)- 
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Figure 5. Schematic illustration of the power and particle balance in a divertor illustrating the roles 
of charge exchange and hydrogen radiation losses from the upstream divertor plasma (away from 
the divertor plate). 
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Figure 6. Ratio of the hydrogen radiated power losses from the upstream divertor plasma to the 
power incident on the divertor plate for ITER and DIII-D conditions. 
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Figure 7. Comparison of Radiation rate coefficients for Be, B, C, Ne and Ar for ADPAK and 
LANL calculations. 
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Figure 8. Radiation loss rate coefficients for Be, B, C, Ne, Ar, Fe, Kr, Mo and W. 
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Figure 9. Normalized Qj^ for the main plasma edge as a function of T e for Be, B, C, Ne, Ar, Fe 
and Kr. 
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Figure 10. Normalized Qn for the divertor plasma and Scrape-Off-Layer for Be, B, C, Ne, Ar and 
Kr as a function of the temperature at the mid-plane separatrix. 
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Figure 1 1 . Normalized Qy for the divertor plasma and Scrape-Off-Layer for C as a function of the 
temperature at the mid-plane separatrix and recycling level (n e x reC y C i e ) 
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Figure 12. Normalized Qy for the divertor plasma and Scrape-Off-Layer for C as a function of the 
temperature at the mid-plane separatrix and neutral fraction no/n e . 
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